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THE DIRECT ALGORITHM FOR SOLVING OF THE GRAPH 
ISOMORPHISM PROBLEM 

RASHIT T. FAIZULLIN, ALEXANDER V. PROLUBNIKOV 



Abstract. We propose an algorithm for solving of the graph isomor- 
phism problem. Also, we introduce the new class of graphs for which 
the graph isomorphism problem can be solved polynomially using the 
algorithm. 



INTRODUCTION 

Complexity of the graph isomorphism problem is not definite still. The graph 
isomorphism problem can not be placed in any of known complexity classes as it 
was stated in pp. It is not proved that the problem is ./VP-complete and there 
is no polynomial-time algorithm for solving of the problem too. Polynomial-time 
algorithms exists only for the special cases of the graph isomorphism problem |3 

Hi. 

We propose an algorithm which is based on the methods of linear algebra. The 
algorithm gives a solution of the problem for graphs from the class that we define 
in this paper. Definition of the class is closely related with graph isomorphism 
notion. We checks weather graphs belongs to the class or not at the process of the 
algorithm operating. We couldn't find a counter example for the algorithm, i.e., we 
couldn't find the pair of isomorphic graphs which can't be defined as isomorphic 
by the algorithm. In particular, at the process of numerous numerical experiments, 
we couldn't find a counter example too. Also, such graphs as regular graphs, which 
traditionally gives hardest cases of the graph isomorphism problem, they belongs 
to this class as numerical experiments shows. 

At the process of the algorithm operating, we implement perturbations of mod- 
ified adjacency matrices of the graphs and solve the systems of linear equations 
associated with them. The algorithm is direct, i.e., the algorithm is not a variation 
of backtracking scheme. There is no search tree, which growth may become uncon- 
trolled at the process of the algorithm operation. The isomorphism of the graphs 
checking at most at n iterations of the algorithm, where n is a number of vertices 
of graphs. If the graphs belong to the class and they are isomorphic, then one 
isomorphism is presented as a result. It is shown that complexity of the presented 
algorithm is equal to 0(n 6 logn) in sense of using elementary operations. 

The paper is organized as follows. In §1, we presents the basic scheme of the 
algorithm and gives a theoretical basis of it using terms of linear algebra. In §2, we 
propose the numerical example of the algorithm implementing. In §3, we consider 
the main computational problems the algorithm deal with; its overall complexity is 
considered. 
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1. THE BASE SCHEME OF THE ALGORITHM 

The graph isomorphism problem. Formulation 1. Suppose Ga = (Va,Ea) 
and Gb = (Vb,Eb) arc two non- weighted non-oricntcd graphs, where Va,Vb arc 
sets of their vertices and Ea,Eb are sets of their edges. \Va\ — \Vb\, \Ea\ = \Eb\- 

Graphs Ga = (Va,Ea) and Gb = (Vb,Eb) are said to be isomorphic if there 
exists a bijection ip : Vb — > Va such that £ Eb {<fi(i) , <f(j)) € Ea- One 

must find this ip or prove that there is no such bijection. □ 

Let us denote isomorphism from graph Ga to Gb as Ga — Gb- Let A be an 
adjacency matrix of graph Ga, i.e., A = (a^), where 



°« \ 



if G 
0, else. 



Let -Bo be an adjacency matrix of Gb- 

We can uniquely set a permutation matrix P v = (py) in correspondence for 
every bijection ip : Vb — * Va- Elements of corresponded permutation matrix P are 
followed: 

1, if i = ip(j), 
0, else. 



Pij 



The following formulation of the graph isomorphism problem is equivalent to the 
first one. 

The graph isomorphism problem. Formulation 2. Suppose Ga = (Va,Ea) 
andGs = (Vb,Eb) are two non-weighted non-oriented graphs. \Va\ = \Vb\, \Ea\ = 
= \Eb\- Let A and B be adjacency matrices of the graphs Ga and Gb- 

One must find a permutation matrix P such that Aq = PBqP^ 1 or prove that 
there is no such permutation matrix. □ 

Let A be an adjacency matrix of Ga = {Va, Ea). Let Da be a diagonal matrix 
of the following form: 

( d& ... \ 
... 



d A 



where 



d£ = J2 a °ij + d = ( h 



and d is a maximal degree of Ga vertices, di is a degree of vertex i £ Va- Matrix 
Db is constructed in similar manner according to matrix Bo. 

The algorithm operates with graph matrices of the following form 

A = A + D Ao , B = B + D Bo - (1) 

A and B are positive definite matrices with diagonal predominance. Let us call 
graph matrices of the form (1) as graph matrices. 

It is clear that the following formulation of the graph isomorphism problem is 
equivalent to the second one and hence it is equivalent to the first one. 
The graph isomorphism problem. Formulation 3. Suppose Ga = (Va,Ea) 
andGs = (Vb, Eb) are two non-weighted non-oriented graphs. \Va\ — \Vb\, \Ea\ = 
= \Eb\- Let A and B be graph matrices of the form (1) for Ga and Gb- 
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One must find a permutation matrix P such that A = PBP~ X or prove that 
there is no such permutation matrix. □ 

Formulation 3 of the graph isomorphism problem may be considered as a re- 
duction of the initial problem to the problem of isomorphism checking of the same 
graphs with additional weighted loops. Matrices of the form (1) are adjacency 
matrices of these graphs. 

Consider the following systems of linear algebraic equations. 

Ax = ej, By = e k , (2) 

where {ej}™ =1 is a standard basis in R™ and A and B are matrices of the form 
(1) (graph matrices) for graphs Ga and Gb- Since A and B are positive definite 
matrices, every system of equations in (2) has solution and the solution is unique. 
Let Xj be a solution of the system of equations Ax = ej and let y k be solution 
of the system of equations By = e k - Vector Xj = (xji, . . . , Xj n ) has the following 
components 

Aij 

where A^j is a cofactor of element of matrix A, i.e. Xj is j-th column of matrix 

A~\ 

Note that A = PBP- 1 iff A' 1 = PB^P' 1 . 

k 

Put A k = A + £j Ej , wh ere 

/ ... ... \ - 1 



E 3 



\ ... ... / - n 

i.e., Ei is a diagonal matrix with the only nonzero element. It is an element at its 
j-th row and it is equal to 1. ej € R, Ej > 0. 

Let G A ' '"' k be the graph which adjacency matrix is matrix A k . 

Put B l = B+ ^2 EjE kj and let Gg' "' 1 be the graph which adjacency matrix is 

j'=i 

matrix B l . Note that B l is defined by sequence {k\, . . . , ki} because we perturb the 
diagonal elements of B by Sj accordingly this sequence while A k is always defined 
by the only sequence {1, . . . , k}. 

Definition 1. We say that graphs Ga — {Va,Ea) and Gb = (Va,Ea) are 
similar and denote it as Ga ~ Gb whenever there exists bijections ip, ipi : Vb — > Va 
such that Xij = y<p(i) Vi {j), = 1, n). □ 

Definition 2. We say that graph Ga is reconstructible whenever 

Ga ~ Gb iff Ga — Gb 

and it holds for every graph Gb- □ 
Let G 1 ^ (i,j = l,n) be the graph which adjacency matrix is Auj\, where A^j^ 
is obtained from matrix A of the form (1) by deleting its i-th row and j-th column. 
If i ^ j, then matrix Au j\ is not symmetrical at general case. G l \ is a weighted 
oriented graph with loops. G\ (i = l,ra) is a graph which is obtained from Ga by 
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deleting vertex i 6 Va and edges which are incident to i. Let us call the graph G^, 
= l,n) as associated graph. 



Let 



k=l k=l 



fc=i /c=i 

where Afc(G^) is fc-th eigenvalue of matrix Auj\, Afc(G^) is fc-th eigenvalue of 
matrix -B(i,j), Afc(Gyi) is fc-th eigenvalue of matrix A, Afc(Gs) is fc-th eigenvalue of 
matrix B. Thus 

lJ detA K A ,Vv detB A B 
This means that graphs Ga and Gs are similar iff there exists bijections ip,<pi ■ 
Vb — > Va such that for graphs G\ and G^ 1 '^' (i,j — l,n) the following holds: 



n A ^ G ^=n A ^ G B^°' ) )- 



There are a lot of ways for constructing cospectral non-isomorphic graphs. Dif- 
ferent types of graph matrices may be considered. But note that there is no uni- 
versal way for this construction. It depends on both the graph representing ma- 
trices and graphs structure. Checking similarity of the graphs, we doesn't check 
their cospectrality. Indeed, we check more slight spectral property of the graphs, 
namely we check the equality of products of eigenvalues from corresponded as- 
sociated graphs spectrums. But we checks coincidence of this characteristics for 
ample quantity of associated graphs at one time. There are n(n — l)/2 of asso- 
ciated graphs for one graph. Checking this characteristics at every iteration of 
the algorithm gives us a powerful method for solving of the graph isomorphism 
problem. We can not find a counter example for the algorithm even though there 
are numerous explorations have been implemented. In particular, there is no non- 
reconstructible graphs in the sense of definition 2 at the library which is located 
at the url http://amalfi.dis.unina.it/graph/ This library is used for testing of the 
most efficient algorithms which are designed for solving of the problem [3]. 

If G A ~ Gb and ip is an isomorphism from Ga to G^, then we have 

/ X n ■ ■ ■ Xlj ■ ■ ■ Xi n \ 



Xnn J 



V Vtp{n)ip(l) ■ ■ ■ y<p(n)ip(j) ■ ■ ■ yip{n)ip{n) / 
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Lemma 1. Let A = PBP^ 1 holds for matrices A and B, where P is a permu- 
tation matrix. If there are no columns of matrix A^ 1 such that they are coincide to 
within the permutation of its components, then P is the only permutation matrix 
for which it holds. 

Proof. The proof of the lemma is obvious. □ 
Remark 1. Clearly, if we take the matrix A itself as a matrix B in the conditions 
of Lemma 1, then we obtain the following statement. If there is no columns of matrix 
A^ 1 such that they are coincide to within the permutation of its components, then 
T(G A ) is trivial (or |r(G A )| = 1). 

Presented algorithm operates with perturbed matrices which are obtained from 
matrices of the form (1). The perturbations implements by adding matrices of the 
form eE^ to the graph matrices. Let A = A, B° = B. At one iteration of the 
algorithm, we consider one matrix A 1 and matrices B J k , where 

A> = A>- 1 + eE\ B 3 k = B^ 1 + eE k . 

Varying k, we finds the vertex kj € Vb and we set kj € Vb in correspondence 
to j e Va iff G A '"' j ~ Gg'"' 3 ~ 1 ' 3 . Perturbations performs until there are no 
columns of A^ 1 that are coinciding within the permutation of its components, i.e., 
perturbations performs until automorphisms groups of the graphs become trivial. 
Let "j <-> kj" denote the correspondence that we set at the process of the algorithm 
implementing. 

The algorithm of graph isomorphism testing 

Base scheme 
Step 0. j := 1, fci := 0. Go to step 1. 

Step 1. If 7 > n, then go to step 7 else chose Sk, I ■— 1, go to step 2. 

Step 2. If I > n, then go to step 6 else if I = k s (1 < s < j — 1), then go to step 4 

else go to step 3. 

Step 3. If C'\ " '' ~ Gg' "' 3 ~ 1 ' , then kj := I and go to step 5 else go to step 4. 
Step 4. I := I + 1. Go to step 2. 
Step 5. j := j + 1. Go to step 1. 

Step 6. Graphs Ga and Gb are not isomorphic or they are not reconstructiblc. 
Stop the algorithm. 

Step 7. Graphs Ga and Gb are reconstructible and isomorphic. Let <p : Vb — > Va 
be a bijection such that f{kj) = j. ip is an isomorphism from Ga to Gb- Stop the 
algorithm. 

Lemma 2. Let P be a permutation matrix and let ip : Vb — * Va be a bijection 
that is corresponded to P. Let A? and B^ be the matrices which are obtained from 
A and B by perturbations of its diagonal elements accordingly to sequences {j}" = i 
and {kj}j- =1 , where ip(kj) = j. 

Then 

A = PBP- 1 Vj A j = PB j p- 1 . 
Proof. A> = Ai- 1 + SjE\ B? = B^ 1 + EjE k K Since 

p{kj)=j^E 3 =PE k ip-\ 

it follows that 

A j = PB 3 P- X & A j - X + SjE j = 
= PiB^ 1 + e j E k >)p- 1 = PB 3 - X P- X + EjPE^p- 1 <^ 
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<S> (A 3 ' 1 = PB } - X P- X and j = <p(kj)). 

The last is true whenever it holds for j = 1, i.e., A = PBP^ 1 . □ 
It follows from lemma 2 that Ga — Gb iff there exists a permutation matrix P 

such that A J = PBiP^ 1 holds for j = l,n whenever fcj = </?(j). I.e., if ~ Gb, 

then, at every j-th iteration, the following must be true: if A J = A 7-1 + E 3 Ei , then 

A j = PB^P- 1 & B j = Bi- 1 +e j E^~ 1 ^. 
By definition, put 

R(A) = {xj}? =1 , R{B) = {y k } n k=1 . 

Definition 3. A set R(A) (R(B)) is called simple if Vj VP ,3/ (/ 7^ j) : x 3 = Px t 
(Vfc VP JM (k 7^ i) yk = Pyi), where P is a permutation matrix. □ 

Note that since graphs matrices have diagonal predominance, the equality Xj = 
Pxi implies that xjj — xu whenever P corresponds to some isomorphism ip. 

At the process of operating of the algorithm, it is necessary to obtain matrices A n 
and B n such that A n = PB n p- 1 and sets R(A n ) and R(B n ) are simple. Thus P 
will be the only permutation matrix that sets the isomorphisms of the graphs that 
are corresponds to the matrices A n and B n . For the graphs that are corresponds 
to the matrices A and B, this matrix will be such permutation matrix too. Recall 
that, for checking the isomorphism of reconstructible graphs, we may check the 
coincidence of vectors Xj and y k for correspondence within to the permutation of 
their components. 

Elements of inverse matrix A -1 are continuous functions of elements of matrix 
A. IiA' = A + ejEi and B' = B + e 3 E k , then 

detA' = det A + ejAjj, dct B' = dct B + e 3 B kk , (3) 

Note that det A', det A, A 3] - , det B', det B, B kk are positive. Again, let P be a per- 
mutation matrix and let ip be a bijection that is corresponded to P. 
Lemma 3. Let 

A = A>- 1 , B = B>- 1 , 

A' = A + eE j , B' k = B + eE k (k = T~R) 

and let Xj be solution of the system of linear equations A'x = e 3 and let y k be 
solution of the system of linear equations B' k y = e k . 

If A = PBP^ 1 holds for some permutation matrix P, then 3\k : x 3 - = Py k . 

Proof. Let x° and yf be solutions of the followed systems of linear equations: 

Ax = ej, By = e t . 

Then 

x1 = ^{A H ,...,A ni ), y? = ^(Bu,...,A nl ). 

Since element a 33 of A and clement b kk of B are the only elements that we change 
by perturbations, then A' ]t = Aj t , {B' k ) kt = B kt (t — 1, n, where {B' k ) kt is a cofactor 
of element b tk of matrix B' k . A and B are matrices with diagonal predominance. 
So if j — (p(k), then x° J3 - — y kk since x° = Py k and ip is an isomorphism from Ga 
to Gb, i.e., A 33 = B kk . Since A = PBP^ 1 , we have detA = det P. Taking into 
account (3), we get det A' = det B' k for every e. Hence if x° } = Py k , then x 3 = Py k 
too because 

_ dct A 1 dot A 

Xj ~ detA 7 X detA™ ' ' ' ' n]) ~ dctA 7 ^"' 
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Vk 



dctB 
dct B' k 



1 



;(Bik, ■ • • , A n k) — 



AetB 
dct B' k 



Vk- 



(4) 



detB 

Let us prove that if columns xj and y k were not coinciding to within the permu- 
tation P of its components before the implemented perturbations, then they won't 
be coinciding after the perturbations too. 

While permutation of matrix rows connected with the same permutation of its 
columns occurs, matrix diagonal elements transforms into transformed matrix di- 
agonal elements. Therefore Xj and yk may coincide to within the permutation of 
its components only if x 33 - = ykk- There are two possibilities of it. First, Ajj = Bkk 
before the implemented perturbations, second, Ajj ^ Bkk before the perturbation. 
At first case, taking into account (3), we get det A' — detB' k . Thus if x° and y k 
arc not coinciding to within the permutation of its components, then Xj and yk are 
not coinciding to within the permutation of its components too since (4) holds. If 
Ajj ^ B kk , then 



_ Ajj 



A 



33 



dct A 1 dct A + eA 



B 



kk 



Bu 



33 



dct B + eBkk 



det B' k 



Vkki 



where e is an arbitrary. Therefore Xj and yk are not coinciding to within the 
permutation of its components. 

□ 

It follows from lemma 3 that the perturbations may be implemented in the way 
such that, in iverse matrix of every graph matrix, there will be no columns that 
are coincides to within the permutation of its components after the perturbations. 
Therefore, it follows from lemma 1 that, at the process of the algorithm operating, 
we transforms the initial graphs into the graphs that have trivial automorphism 
group. Similarity of the graphs is preserving while perturbations implements. In 
case of reconstructible graphs, it follows from lemma 2 that correspondence j <-> kj 
is an isomorphism, where j <-> kj is a correspondence that is settled at the iterations 
of the algorithm. 



2. EXAMPLE OF THE ALGORITHM IMPLEMENTING 

Let us consider example of the algorithm implementing. Let Ga and Gb be 
isomorphic tested graphs. For these graphs, matrices of the form (1) are followed 
matrices: 



/ 8 1 1 1 1 \ 

1 7 1 1 

1 1 7 1 

1 7 1 1 

10 17 1 

\ 1 1 1 1 8 / 



B 



( 8 1 1 1 1 \ 

8 1111 
117 10 
1117 

1 1 7 1 
\1 10017/ 
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Let machine numbers mantissas length is equals to the number of graph vertices 
(n = 6). Then 



(AT 1 



/ 8 
1 
1 
1 
1 

V o 

0.134 
-0.018 
-0.018 
-0.018 
-0.018 
\ 0.008929 



1111 
7 10 
17 
7 1 
17 
1111 



\ 
1 
1 
1 
1 

8 / 



-0.018 

0.15 
-0.016 



-0.018 
-0.016 
0.15 



-0.018 -0.018 0.008929 \ 
0.004464 0.004464 -0.018 
0.004464 0.004464 



0.004464 0.004464 
0.004464 0.004464 
-0.018 -0.018 



0.15 
-0.016 
-0.018 



-0.016 
0.15 
-0.018 



-0.018 
-0.018 
-0.018 
0.134 



/ 



and 



o\-i - 







1 \ 



/ 8 

1 
1 
1 

V 1 

0.134 
0.008929 
-0.018 
-0.018 
-0.018 
\ -0.018 

In this case, possible solutions of the graph isomorphism problem are permuta- 
tion matrices which are corresponded to the following substitutions: 



/ 



1 1 
111 
1 
7 
1 7 1 
10 17/ 

0.008929 - 
0.134 
-0.018 
-0.018 
-0.018 
-0.018 



-0.018 
-0.018 

0.15 
-0.016 



-0.018 
-0.018 
-0.016 
0.15 



-0.018 
-0.018 



-0.018 \ 
-0.018 



0.004464 0.004464 
0.004464 0.004464 



0.004464 0.004464 

0.004464 0.004464 

0.15 -0.016 
-0.016 0.15 



1 


2 


3 


4 


5 


6 


1 


3 


4 


5 


6 


2 


1 


2 


3 


4 


5 


6 


1 


4 


3 


6 


5 


2 


1 


2 


3 


4 


5 


6 


1 


6 


5 


3 


4 


2 


1 


2 


3 


4 


5 


6 


2 


3 


4 


6 


5 


1 


1 


2 


3 


4 


5 


6 


2 


5 


6 


3 


4 


1 



Let 



P' 



1 2 


3 


4 


5 


6 


1 3 


4 


6 


5 


2 


1 2 


3 


4 


5 


6 


1 5 


6 


3 


4 


2 


1 2 


3 


4 


5 


6 


1 6 


5 


4 


3 


2 


1 2 


3 


4 


5 


6 


2 4 


3 


5 


6 


1 


1 2 


3 


4 


5 


6 


2 5 


6 


4 


3 


1 


1 2 


3 


4 


5 


6 


2 6 


5 


4 


3 


1 




2 


3 


4 


5 


(I 


4 


5 


3 


6 



1 


2 


3 


4 


5 


6 


1 


4 


3 


5 


6 


2 


1 


2 


3 


4 


5 


6 


1 


5 


6 


4 


3 


2 


1 


2 


3 


4 


5 


6 


2 


3 


4 


5 


6 


1 


1 


2 


3 


4 


5 


6 


2 


4 


3 


6 


5 


1 


1 


2 


3 


4 


5 


6 


2 


6 


5 


3 


4 


1 
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Note that, in particular, first column of B^ 1 coincide to within the permutation 
P 1 with first column of A^ 1 but P' is not sets isomorphism from graph Ga to 
graph Gb- One can see that P' won't be such matrix after the 3-rd iteration of the 
algorithm. 

Iteration 1. Let e\ = 0.1; then 



( i 


i + £l 


1 


1 


1 


1 


o \ 




1 


7 


1 








1 




1 


1 


7 








1 




1 








7 


1 


1 




1 








1 


7 


1 


{ 





1 


1 


1 


1 


8 J 



0.132 
-0.018 
-0.018 
-0.018 
-0.018 
0.008811 



-0.018 
0.15 

-0.016 
0.004433 
0.004433 

-0.018 



-0.018 
-0.016 

0.15 
0.004433 
0.004433 
-0.018 



-0.018 
0.004433 
0.004433 

0.15 
-0.016 
-0.018 



-0.018 
0.004433 
0.004433 

-0.016 
0.15 

-0.018 



0.008811 
-0.018 
-0.018 
-0.018 
-0.018 
0.134 



thus if we put k\ = 1, then we obtain 



( 8 


+ ei 





1 


1 


1 


1 \ 







8 


1 


1 


1 


1 




1 


1 


7 


1 










1 


1 


1 


7 










1 


1 








7 


1 


V 


1 


1 








1 


7/ 



0.132 
0.008811 
-0.018 
-0.018 
-0.018 
-0.018 



0.008811 
0.134 
-0.018 
-0.018 
-0.018 
-0.018 



-0.018 
-0.018 

0.15 
-0.016 
0.004433 
0.004433 



-0.018 
-0.018 
-0.016 

0.15 
0.004433 
0.004433 



-0.018 
-0.018 
0.004433 
0.004433 

0.15 
-0.016 



-0.018 
-0.018 
0.004433 
0.004433 
-0.016 
0.15 



Correspondence j <-> kj has the following form after iteration 1: 



1 2 3 4 5 6 
1 



where points denote correspondences which are not settled yet. 
Iteration 2. Let us put e-i = 0.2 
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( 8 


+ ei 


1 


1 


1 


1 


o \ 




1 


7 + £ 2 


1 








1 




1 


1 


7 








1 




1 








7 


1 


1 




1 








1 


7 


1 


V 





1 


1 


1 


1 


8 / 



0.134 


-0.017 


-0.018 


-0.018 


-0.018 


0.00875 


-0.017 


0.146 


-0.016 


0.004303 


0.004303 


-0.017 


-0.018 


-0.016 


0.15 


0.004447 


0.004447 


-0.018 


-0.018 


0.004303 


0.004447 


0.15 


-0.016 


-0.018 


-0.018 


0.004303 


0.004447 


-0.016 


0.15 


-0.018 


0.00875 


-0.017 


-0.018 


-0.018 


-0.018 


0.134 



(B ) = 



( 


3 + £i 





1 


1 


1 


1 \ 







8 


1 


1 


1 


1 




1 


1 


7 + £ 2 


1 










1 


1 


1 


7 










1 


1 








7 


1 


V 


1 


1 








1 


7/ 



0.132 


0.00875 


-0.017 


-0.018 


-0.018 


-0.018 


0.00875 


0.134 


-0.017 


-0.018 


-0.018 


-0.018 


-0.017 


-0.017 


0.146 


-0.016 


0.004303 


0.004303 


-0.018 


-0.018 


-0.016 


0.15 


0.004447 


0.004447 


-0.018 


-0.018 


0.004303 


0.004447 


0.15 


-0.016 


-0.018 


-0.018 


0.004303 


0.004447 


-0.016 


0.15 



Correspondence j <-> kj has the following form after iteration 2: 



1 2 3 4 5 6 
13.... 



Iteration 3. £ 3 = 0.3 



V 



f £l 


1 


1 


1 


1 





1 


7 + £ 2 


1 








1 


1 


1 


7 + £3 








1 


1 








7 


1 


1 


1 








1 


7 


1 





1 


1 


1 


1 


8 



/ 



/ 0.132 
-0.017 
-0.017 
-0.018 
-0.018 

\ 0.008659 



-0.017 
0.146 

-0.015 
0.004324 
0.004324 

-0.017 



-0.017 
-0.015 
0.144 
0.004255 
0.004255 
-0.017 



-0.018 
0.004324 
0.004255 

0.15 
-0.016 
-0.018 



-0.018 
0.004324 
0.004255 

-0.016 
0.15 

-0.018 



0.008659 \ 
-0.017 
-0.017 
-0.018 
-0.018 
0.134 / 
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(B3)-l = 



/ 1 


i + £l 





1 


1 


1 


1 \ 







8 


1 


1 


1 


1 




1 


1 


7 + £ 2 


1 










1 


1 


1 


7 + £.3 










1 


1 








7 


1 




1 


1 








1 


7 / 



0.132 
0.008659 
-0.017 
-0.017 
-0.018 
-0.018 



0.008659 
0.134 
-0.017 
-0.017 
-0.018 
-0.018 



-0.017 
-0.017 
0.146 
-0.015 
0.004324 
0.004324 



-0.017 
-0.017 
-0.015 
0.144 
0.004255 
0.004255 



-0.018 
-0.018 
0.004324 
0.004255 

0.15 
-0.016 



-0.018 
-0.018 
0.004324 
0.004255 
-0.016 
0.15 



Correspondence j <-» kj has the following form after iteration 3: 



1 2 3 4 5 6 
13 4. . . 



Iteration 4. £4 = 0.4 



( 0.132 -0.017 
-0.017 0.146 

-0.017 -0.015 

-0.017 0.004079 

-0.018 0.004351 

\ 0.008541 -0.017 



/ 


3 + £l 


1 


1 


1 


1 


\ 




1 


7 + £ 2 


1 








1 




1 


1 


7 + £3 








1 




1 








7 + e 4 


1 


1 




1 








1 


7 


1 


V 





1 


1 


1 


1 


8 / 



/ 8 + £i 



V 



-0.017 
-0.015 
0.144 
0.004014 
0.004282 
-0.017 

1 



8 1 1 

1 7 + £ 2 1 

1 1 7 + £ 3 
1 
1 



0.132 
0.008541 
-0.017 
-0.017 
-0.017 
-0.018 



0.008541 
0.134 
-0.017 
-0.017 
-0.017 
-0.018 



-0.017 
-0.017 
0.146 
-0.015 
0.004079 
0.004351 



-0.017 
0.004079 
0.004014 
0.142 
-0.015 
-0.017 

1 
1 




7 + £ 4 
1 

-0.017 
-0.017 
-0.015 
0.144 
0.004014 
0.004282 



-0.018 
0.004351 
0.004282 

-0.015 
0.15 

-0.018 

1 \ 
1 




1 

-0.017 
-0.017 
0.004079 
0.004014 
0.142 
-0.015 



0.008929 \ 
-0.017 
-0.017 
-0.017 
-0.018 
0.134 / 



-0.018 
-0.018 
0.004351 
0.004282 
-0.015 
0.15 



Correspondence j <-> kj has the following form after iteration 4: 

1 2 3 4 5 6 
1 3 4 5. . 
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The sets R{A 4 ) and R(B A ) are already simple before 5-th iteration. Thus we 
may not perturb matrices A 4 and B 4 no more and we can set the correspondence 
for the rest of the graphs vertices (j = 5, 6). 

As a result, we obtain the isomorphism from Ga to Gg. It is 

/ 1 2 3 4 5 6 \ 
^ 1 3 4 5 6 2 ) ' 

3. COMPUTATIONAL EFFICIENCY OF THE ALGORITHM 

3.1. Localization of separated solutions of linear equations systems. Let 

us denote solution of system A>x = as x^\ By definition, put 

UA r ) = {xV> \Vp,q:l<p,q<k3P: = Px&}. 

We have R(A) = URi(A). Let m(A j ) be the amount of sets Ri{A j ) for matrix A j . 
R(A) is simple iff |i?j(A)| = 1 for every i and hence m(A) = n. Let us consider 
these sets at the process of the perturbations of matrix A. 

Definition 4. Let us say that set Ri(A^ r ) is splitted by separation of some 
vector Xj* -1 ^ G i£i(A J_1 ) from it at j-th iteration of the algorithm if x^p G i?fc(A J ) 
for some fc and |iifc(.A J ')| = 1. Let us call such vector as a separated vector. □ 

If we get m(A : >) > m{A^~ 1 ) and, finally, m(A n °) = m(B n °) = n for some no 
(no < n) while perturbing the graphs matrices at iterations of the algorithm, then 
sets R(A) and R(B) will be simple sets R{A n °) and R(B n °) at some no-th iteration. 

By definition, put 

4? = 141 -zgij. 4? e i? fe (^), 4f g (5) 

Let us understand A$ as a distance between different sets Rk (A^ ) and Ri (A* ) . If 

^ x^, then 4f G #fc(-4) and 4f G fl;(A), where fc ^ I. 
By definition, put 

4? = l41 - 4?U 4? e J2 fc (^), 4f G Rt(A*). (5) 

may be understand as a distance between different sets Rk(A^) and Ri(Ai). 

If 4f fel ^ 4iL' then x fe? e ^O 4 ) and 4f e R ii A ), where fc ^ /. 

At the presented above example, the splitting of sets may be illustrated as follows. 
Before 1-st iteration of the algorithm: 

Ri (a°) - {4°\4 0) }, m*°) = {4 0) ,4°\4 0) ,4 0) }, 

= 0.016, 

m(A°) = 2. 

After 1-st iteration: 

RiiA 1 ) = {x[ 1} }, rm 1 ) = {4 1) ,4 1) 5 4 1) .4 1) }. ma 1 ) = {4 1} i: 

A$ = 0.018, A$ = 0.002, A$ = 0.016, 
miA 1 ) = 3. 

After 2-nd iteration: 

R,{A 2 ) = {4 2 >}> MA 2 ) = {4 2) >, MA 2 ) = {x {2) }, 
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R 4 (A 2 ) = {xf\xf ) }, R 5 (A 2 ) = {x^}, 
A$ = 0.014, A$ = 0.016, A^ = 0.018, A$ = 0.002, A$ = 0.004, 
A^ } = 0.004A^ } = 0.012, A$ = 0, A^ = 0.016, A 4 2 5 } = 0.016, 

m(A 2 ) = 5. 

After 3-rd iteration: 

RM 3 ) = {4 3) }, MA 3 ) = {x {3) }, MA 3 ) = {x {3) }, 
R,(A 3 ) = {x (3 \x i3) }, R 5 (A 3 )={x i3) }, 
A^ } = 0.014, = 0.012, A^ = 0.018, A^ = 0.002, A^ = 0.002, 
A^ } = 0.004, A^ = 0.012, A^ = 0.006, A^ } = 0.01, A^f = 0.016, 

m(A 3 ) = 5. 

After 4-th iteration: 

R,{A 4 ) = {4 4) L MA 4 ) = {x (4) }, R 3 (A i ) = {4 4) }, 

R 4 (A 4 ) = {4 4) l, MA 4 ) = {x {4) }, RM 4 ) = {4 4) l- 
A$ = 0.014, A$ - 0.012, A^ } = 0.001, A$ = 0.018, A$ = 0.002, 
A$ = 0.002, A$ = 0.004, A^ = 0.004, A$ = 0.012, A^ } = 0.008, 
A^ } = 0.006, A$ = 0.01, A 4 5 } = 0.008, A$ = 0.008, A^ = 0.016, 

m(A 4 ) = 6. 

Finally, we obtain that every initial set of vectors Ri{A) is splitted and uq = 4 
whereas n — 6. Note that there are no separated vectors after 3-rd iteration but all 
of column- vectors Xj (J = l,n) are separated from each other after 4-th iteration. 
Matrix conditional number is the following value: 

• ||Ac||||g|| 

c / ;i#oiP£iiiMi 

where || • || is Euclidian norm of vectors in R™. Let Amax be a maximal eigenvalue 
of matrix A spectrum and let A mm be a minimal eigenvalue of matrix A spectrum. 

sup ||Ac||/||^|| 
infPeil/Heil A min 

A mm =/= because A is a positive definite matrix. 

Let us consider the system of equations (A + C)y = / that is obtained from 
Ax = f by perturbation of matrix A by matrix C . Let 

\\C\\ <f) 

FiT 

where 

n ll^ll 

\\A\\ = Amax for positive definite matrix A The following theorem holds UJ. 
Theorem (Godunov). If 9(J,(A) < 1, then 

\\y-x\\ < en{A) 



n{A) — sup 



1 - 0n(A) ' 
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□ 

The following holds for symmetrical matrix conditional number 0: 

ri(A) 



x(AY 



where r)(A) = max (an + x(A) = min (an - \ a ij\)- Since a kk = 

l<i<n jjti l<t<n 

= d + d k + £fe at j-th iteration, where is value of perturbation of fc-th matrix A 
diagonal element and £fc = if k < j. Let i\ be number of the row that gives rj(A) 
and let ii be number of the row that gives x(A). We have 

n 

■q(A) = a ilh + loijjl = d + e hll +d n +d h = d + e hll +d + d= 3d + e hil , 



= a l2l2 + E i3 i 2 - 7 ^ |ai 2 j| = d + e iin + d i2 - d l2 = d 

Therefore 



,(A)<^- M - 



X(A) d + s i2i2 ' 
If Ejjij = l/n Pl , £i 2 i 2 = l/n P2 , where pi, p2 £ N, then 

U) < gd+fi^ = 3d+l/^ = nP2 _ pi < nP2 _ pi Adn^ = ^ 

v d + e i2 i 2 d+l/nP? dnP^ + l ~ dn? 2 

So the stated above theorem take the following form in our case: 
If 9 < 1/4, then 

lly-sHflj^NI. (6) 

□ 

The perturbations of matrices must be small enough for guaranteed isolation of 
vector Xj* separated from Ri(A^ r ) at j-th iteration. This means that if Xj ^ is 

separating from Ri(A^~ 1 ^), then it is necessary that, after the perturbation, X ' 
will be sufficiently far from other sets in sense of introduced distance A (5). We 
can guarantee this if there is a localization of the linear equations systems solutions 
take place. Here, we means localization about the splitted set iZ,-(AW — 1 ') in the 
sense of distance A. 

It follows from (6) that 

Vz: |U«_ s (*-i)||< «C^|uy-i)|i 



1 - 40( £j 



where 

l|C|| _ \\ Ej E*\ 



WAi^W WAi-m \\Ai- 

So we can put 



In- 



putting £j = l/n p , we get 



H_(j-i) _ t «)ii < A / nP || r tf-i)|| - 4 'V nP i| r 0-i)i 
11 * 1 11 \\AJ\\-4/nP ]l < 11 pJ||-4.1/np" ' 1 
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We have 

^{A j - r ) = A m ax/A min < 4, 

where Amax is a maximal eigenvalue of matrix spectrum and A m j n is a mini- 
mal eigenvalue of matrix A 3 spectrum. Since Gershgorin's Theorem 0, we have 

Amax > d. 

Taking into account that A m j n < we get 

\\A j - x \\ > A min > A m ax/4 > d/4. 

So we have 

ll^-^IK^- 13 !!- 

Since A m j n > 1, ||ej|| = 1 and | = A~ 1 ej, then < 1. So we have 

II^IK^- 



Therefore it follows from (7) that 



Vi: \\x^-xf\\< 



(PnP- 1 ' 
It is easy to prove that 



So we have 



V- ^-^K^I (8) 

(8) gives an estimate of localization of vectors from sets Ri(A^) at the process of 
the algorithm iteration. 

3.2. Splitting of the sets of linear equations systems solutions. Imple- 
mented perturbations must be sufficient for splitting of the sets Rk at the process of 
operation of the algorithm. The algorithm is computationally effective only if, for 
every j, fc, I, distances A.$ will be a non-zero machine number. Let us show that 
we can obtain adequate accuracy using machine numbers which mantissas length 
not exceeding n. 

Let x'f 1 , Xj' € Ri(A^ 1 ) for some i, i.e., there are at least two elements 
that prevents us to set the univocal correspondence between the vertices of graphs 
at j-th iteration. If ^ Xj) after j-th iteration, then there is no permutation 

matrix P such that x[^ = Px^ . Let us estimate the lower bound of the distance 

(i) (i) 

between such x{ and Xj ' after j-th iteration of the algorithm. 

Theorem 1. If there exists permutation matrix P such that x^ = Pxf ^ 
and < Sj < 1, then 

aM _ \ Jj) _ Jj) \ > \ 
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Proof. Put A = Ai- 1 , A' = A>, x t = x\ = xf \ Xj = xf 1} , = xf. 

If x'a — Px' it then x'- • = x' u . 



A'- ■ A' 

31 x { _ A H 



" det A" M detA r 
We have 

det A' del .1 • £j -Ajj , A^ = Ajj 

because the diagonal element cijj is the only element of A that changes at the 
iteration. Since 



we get 



A' A' 



A, ;i — An + EjAnjj, 



A ' ■ .1 . i 

33 ^Hi T fc J-^MJJ 



det A' det A' 



det A' det A' 



det A + EjAjj 



Anjj is a determinant of the submatrix of A that is obtained by deleting of 
i-th row and i-th column and j-th row and j-th column. A is a symmetrical 
positive definite matrix with a diagonal predominance hence the following holds for 
A (Hadamard): 

H k = \dkk\ - X! I afei l = dk > °' i = 1 > ri ' 

where 

, = f d + e fc , if 1 < k < j, 
k \ d, if j < k < n. 

Consequently (see |S]), 

Aajj > Hi ...Hi.. .Hj .. .H n = d[. . .dj. . .d] .. .d' n > d" (10) 

where Hi means that Hi is absent at the product. 

On the other hand, there is a following estimate for the maximal eigenvalue 
Amax of A (Gershgorin's Theorem, see [7]): 

Amax < 3d + max e k < 3d + 1, 

i<fe<i 

if Vfc £fc < 1 . Therefore 

71 

det A = Y[ X k < A max < (3d + 1)", (11) 

k=l 

n-1 

Ayj = I] X 'k < ASTi < (3d + I)"" 1 , (12) 
fc=i 

Let A max be the maximal eigenvalue of the matrix that is remaining from the 
matrix A by deleting its i-th and j-th rows and columns. A max — Amax E!> 
consequently, combining (10), (11), (12), from (9) we get: 

\rJ _ T ' I _ £ 3^ii,33 > E 3& n £]d n 

1 " iil detA + ^Ajj {3d+l) n + e J (3d+l) n ~ 1 {Sd)^ 1 +e,{3d) nl 
i.e., 



A ? " l^i? ^ I > 3 n d 2 {M/£j + 1} ■ ( 13 ) 

□ 
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Let we shall be solve the systems of linear equations A>x = eu and B^y = eu 
(j, k = 1, n) by Gauss-Seidel iterative method. Actually, it is doesn't really matter 
what method of solution is chosen: we just need a proper accuracy of solutions 
which we can obtain in a proper time. But Gauss-Seidel method also gives us an 
efficient way to estimate the needed machine number mantissas length. 

The following theorem may be found in PJ. 

Theorem 2. Let ^ < 7|<2«|)7 < L for every i = l,n. Then 

||xi-£|||<7l|xi-5r 1 H. 
where Xi is exact solution of the system of linear equations A* x — e^, if is an 
approximate solution of the system of equations at s-th iteration of Gauss-Seidel 
method. □ 
If A = A? for some j, then 7 < 1/2. Therefore at s-th iteration of Gauss-Seidel 
method we have: 

\x u - x s u \ < \\ Xi - xt\\ < —6°, \xjj - x] 3 \ < \\xj - x^W < —5°, 

where Sq is a mistake of initial approximation. If Ej = \/n p at j-th iteration, then, 
taking into account (12), we get 

I CO _ CO I > i 

I a a I ^ 3n d 2( MnP + 1 y 

Hence if 

— 5° < - (14) 

2 s-i - 2, n d 2 {MnP + iy 1 ' 

then difference \Xjj — x^ \ may be fixed by machine numbers which mantissas 
length is restricted by n. Let us estimate the number of Gauss-Seidel method 
iterations s that must be implemented to fix this value. Using (14), we get 

3 n d 2 {3dn p + 1)5° <2 S ~ 1 . (15) 

Taking the logarithm of (15), we obtain 

nlog 2 3 + log 2 d 2 {3dn p + 1) + log 2 5° < s - 1, 
i.e., the number of needed method's iterations may be estimated as 

s > n log 2 3 + p log 2 n + 3 log 2 d + log 2 5 + 3. (16) 

Let p > be an integer that defines the value of matrix diagonal clement per- 
turbation: £j — l/n p . Let Xj ^ be a vector that isolates from set R^A^ 1 ) at 
j-th iteration. 

Let us choose p > such that 

4 A, 
< — , 



d 2 n p 



i.e., 



p>1 ° Sn d 2 ~A~ + 1, (17) 



where A,- = min A$ . It follows from (7) that 



k^l 



Vz: |4- 1) -4 ) |<-^<^ 
1 " d 2 n p n 
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And, using of theorem 1, we have 



Vi: - m -t<|^ ) -4 ) |< 



3 n d 2 (3d/e J - + 1) 1 " " 1 d 2 n?> 

after j-th iteration. This means that, at j-th iteration, we can choose such per- 
turbation that will give us a separation that is fixable by machine numbers with 
restricted mantissas length and localization of vectors Xi at the same time. This 
mean that we can set the interval in which the separating of vectors occurs. We 
can perform this process in the way such that if we isolating some vector xj from 
then this vector won't be placed in another such set of vectors at the 
next iteration unless power of this set is equal to 1. 
Lemma 4. 

1 

j > 3 n d 2 (3dn n + 1) 

and 

Sj < (n + 3) log 2 n 

at every j-th iteration of the algorithm. 

Proof. Let us consider the following case of the graph isomorphism problem. Let 
graphs Ga and Gb be complete graphs, i.e., every vertex of each graph is adjacent 
to every other vertex of the graph. It is clear that checking of the isomorphism of 
complete graphs is easy taken by itself since every complete graph is isomorphic 
to another complete graph on the same vertices set. The purpose of introducing 
the following example is to show the implementation of the algorithm in hardest 
case because m(A) = 1, |i?i(A)| = n and d — n here, i.e d is maximal for every 
non- weighted graph, so, taking into account theorem 1, Aj take on a minimal value 
in this case. 

1. Put e\ = 1/n at 1-st iteration, x^ is separated at this iteration and, for every 

l/(3 n d 2 (3rfn+ 1)) < |ccn -XjA < 4/(d 2 n). 
Needed number of Gauss-Seidel method iterations is s, where 
s = n log 2 3 + log 2 n + 3 log 2 d + log 2 5° + 3. 

2. Put e 2 = 1/n 2 at 2-nd iteration. 2-nd vector is separated at this iteration and, 
for every j ^ 2, 

l/(3 n d 2 (3dn 2 + 1)) < |x 22 - Xjj \ < 4/(d 2 n 2 ). 

Needed number of Gauss-Seidel method iterations is s, where 

s = n log 2 3 + 2 log 2 n + 3 log 2 d + log 2 5° + 3. 

Continuing separating solutions of linear equations systems in this way, we obtain 
£2 = 1/n" -1 at (n — l)-th iteration of the algorithm and (n — l)-th vector is 
separated and, for every j / n - 1, 

l/(3"d 2 (3dn"- 1 + 1)) < |x„_i, n _i - Xjj\ < V^n"" 1 ). 

Needed number of Gauss-Seidel method iterations is s, where 

s = n log 2 3 + (n — 1) log 2 n + 3 log 2 d + log 2 6 + 3. 

Since exact solution of the systems belongs to segment [0, 1]™ G M™, we may put 
5° = 0(y/n) and, consequently, s < (n + 3) log 2 n. □ 
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As a result, if we set machine numbers mantissas length is equal to n and perform 
s = (n + 3) log 2 n Gauss-Seidel method iterations for solution of every system of 
linear equations at iterations of the algorithm, then we obtain that m(A n ^ 1 ) = n 
while the algorithm is computationally effective to fix occurring split of sets of linear 
equations systems solutions. If we use the extended numerical machine numbers 
type that realized in languages such as Object Pascal and C++, then we can 
operate with machine numbers at the range from 3, 6 x io~ 4951 to 1,1 x 10 4932 . 
Thus the algorithm may be efficient for reconstructible graphs for which the number 
of vertices is lesser than 4951 and it is possible to operate with standard extended 
machine numbers without any additional procedures for correct operating with 
them. 

3.3. Overall complexity of the algorithm. At every j-th iteration of the algo- 
rithm, we checks similarity of G 1 ^ and Gg'"' at step 3, where I can take 
on a value from 1 to n. Checking of similarity implements using inverse matrices 
to matrices A 7 and B l (I = l,n), where 

j-i 

i=l 

We obtain inverse matrices solving the systems of linear equations. Thus we must 
solve n+nxn linear equations systems at most. Let N$ be the number of elementary 
machine operations that we must implement to solve linear equations system with 
needed accuracy. Then complexity of finding all needed inverse matrices is equal 
to n(n + 1)N S . 

Checking of coincidence to within the permutations of inverse matrices columns 
can be implemented at 0(n log n) elementary machine operations. To check the 
similarity of two graphs G 1 ^ " J and Gg''"' we must implement checking of n 2 
pairs of inverse matrices columns at most. Therefore complexity of checking for 
similarity of two graphs is equal to 0(n 3 logn). The number of such graph pairs 
may be equal to n at most. 

Hence overall complexity of one iteration of the algorithm is equal to 

0(n(n + 1)N S ) + n ■ 0(n 3 logn) = 0(n 2 N s + n 2 logn). (18) 

Let us estimate the value of N$ for Gauss-Seidel method of solution of linear 
algebraic equations. 

At j-th iteration of the algorithm, we calculate p such that 4/(<i 2 n p ) < Aj/n, 
where Aj — minA^, i.e., 

p=iog "rf4- +i - 

Then we implement s iterations of Gauss-Seidel method, where, from lemma 4, we 
have 

s = (n + 3) log 2 n. 

Let Ngj be a complexity linear equations systems solution at j-th iteration. Then 

N 3 S = 0(n 2 ) ■ (n + 3) log 2 n 

since one iteration of Gauss-Seidel method has complexity that is equal to 0{n 2 ). 
Thus, using (18), summary complexity of all iterations of the algorithm may be 
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estimated as 



0(n 2 ) • 0{n 2 ) ■ (n + 3) log 2 n = 0(n 6 + n l 




(19) 



Complexity of the algorithm may be reduced if the sets R(A n °) and R(B n °) 
become simple at no-th iteration, where n$ is much lesser than n. In particular, 
applying the algorithm for isomorphism checking of regular graphs which are lat- 
tices on torus, the sets R(A n °) and R(B n °) become simple when no = *Jn. Vertex 
degree is equal to 4 at this case. Let us remark that regular graph isomorphism 
testing is one of the hardest case of the graph isomorphism problem for most ef- 
ficient algorithms which are designed for solving of the problem. As an example, 
complexity of NAUTY algorithm JU| becomes exponential when vertex degree of 
regular graphs is equal to 4 [TT], 

Remark 2. Since (11) and (12) are obtained using majorization of graph matrix 
A eigenvalues by A m ax, then (13) holds regardless of multiplicity of form (1) matrix 
eigenvalues and (13) holds regardless of multiplicity of graph adjacency matrix 
eigenvalues too. Thus (20) holds regardless of multiplicity of of reconstructible 
graphs adjacency matrix eigenvalues. 

Remark 3. Since d < n, we have 



Thus (20) holds regardless of maximal degree of reconstructible graphs vertices. 



An algorithm for solving of the graph isomorphism problem is presented. The 
algorithm solve the problem for the class of reconstructible graphs which is de- 
fines at the paper. The algorithm checks weather graphs belongs to the class of 
reconstructible graphs or not at the process of algorithm implementing. 

Proved the theorem that gives an estimate of computational complexity of the 
algorithm. It is shown that algorithm is polynomial in sense of using elementary 
machine operations and in the sense of used memory too since used machine num- 
bers mantissas length is restricted by polynomial of number of graph vertices. It is 
shown that this is holds regardless of maximal degree of the graphs, graphs genus, 
graph eigenvalue multiplicity etc., i.e., using the algorithm, solution of the graph 
isomorphism problem has no specific that may be determinated by any graph char- 
acteristics that is usually considered. 
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